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Using benzene sandwiched between two Au leads as a model system, we investigate from first 
principles the change in molecular conductance caused by different atomic structures around the 
metal-molecule contact. Our motivation is the variable situations that may arise in break junction 
experiments; our approach is a combined density functional theory and Green function technique. 
We focus on effects caused by (1) the presence of an additional Au atom at the contact and (2) 
possible changes in the molecule-lead separation. The effects of contact atomic relaxation and two 
different lead orientations are fully considered. We find that the presence of an additional Au atom at 
each of the two contacts will increase the equilibrium conductance by up to two orders of magnitude 
regardless of either the lead orientation or different group- VI anchoring atoms. This is due to a 
LUMO-like resonance peak near the Fermi energy. In the non-equilibrium properties, the resonance 
peak manifests itself in a large negative differential conductance. We find that the dependence of 
the equilibrium conductance on the molecule-lead separation can be quite subtle: either very weak 
or very strong depending on the separation regime. 

PACS numbers: 73.40.Cg, 72.10.-d, 85.65.+h 



I. INTRODUCTION 

Understanding electron transport through nanoscale 
junctions or molecular devices connected to metallic elec- 
trodes may be the basis of future molecular electronics 
technology. 1,2,3,4,5 One of the critical issues in this regard 
is to construct contact structures which can provide both 
useful stability and high contact transparency. 

In many recent experiments^SiL& Au electrodes were 
used as leads for transport measurements because of the 
high conductivity, stability, and well-defined fabrication 
techniques involved. A common way to construct a lead- 
molecule-lead (LML) system is by using a mechanically 
controllable break junction (MCB). These can be made 
either through direct mechanical meansSiSiS or through 
electromigrationifiiiiiiSii^ In these break-junction exper- 
iments the detailed atomic structure of the molecule-lead 
contacts of a LML system is unknown. In fact, because of 
the atomic scale roughness of the break surface, different 
atomic scale structures of the contact may occur in dif- 
ferent experiments. Up to now, it has been difficult to in- 
vestigate and control experimentally the detailed contact 
atomic structure and find out its influence on the electron 
transport through a molecular device. Here theoretical 
modeling/simulation free from empirical parameters can 
play an important role in understanding, interpreting ob- 
served experimental behaviors, or doing pre-designs for 
good contact structures. 

Because of the need for an atomic scale descrip- 
tion, ab initio density functional theory (DFT) 1 ^ 5 - 
is a natural approach to molecular electronics. One 
route works explicitly with scattering states via the 
Lippmann-Schwinger equation^SiiLi2iiS The more 
common method, however, is to combine DFT using 
a localized basis set for molecular electronic struc- 
ture with the non-equilibrium Green function method 



(NEGF) 2 ^ for electron transport. 22 i 23 i 24 i 25 i 26 i 27 In 
this method some parts of the electrodes of a LML 
system can be included in the device region to form 
an "extended molecule", and therefore the specific 
contact atomic structure and relaxation can be fully 
considered in principle although in practice the atomic 
structure is usually predetermined to avoid the heavy 
computational effort. 28 i 29 i 30 i 31 i 32 i 33 i 34 i 35 i 36 i 37 i 38 i 39 i 40 
In the implementations of this approach, some 
researcher a 22 i 23 i 24 i 28 i 29 i 30 i 31 i 32 adopted quantum chem- 
istry methods for the DFT calculation in which a 
cluster geometry is used for the device region, while 
othera 25 i 26 i 27 i 33 i 34 i 35 i 36 i 37 i 38 i 39 i 40 used a periodic geome- 
try (as in solid state physics) for the device region. The 
advantages of the latter are that the electronic structure 
of the device region and the two leads can be easily 
treated on the same footing and the infinite LML system 
is nearly perfect in geometry without any artificially 
introduced surface effect. 

Previously we developed a self-consistent approach 
within the DFT+NEGF method for calculating electron 
transport through molecular devices. 27 Our approach is 
simple while strict: the non-equilibrium condition un- 
der a bias is fully included in the NEGF rather than the 
DFT part. Therefore, it is straightforward to combine 
with any electronic structure method that uses a local- 
ized basis set. More importantly, in this way the problem 
of solving for the Hartree potential under a bias field with 
unphysical potential jumps at the two boundaries of the 
DFT supercell is avoided. In our method large parts of 
the two metallic leads of a system are included in the 
device region so that the molecule-lead interactions (in- 
cluding electron transfer and atomic relaxation) are fully 
included, and the electronic structures of the molecule 
and the two leads are treated exactly on the same foot- 
ing. 
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FIG. 1: Optimized atomic structures of the LML systems investigated, obtained by relaxing fully the molecule, the first two 
atomic layers of the two lead surfaces, and the molecule- lead separation, (a), (b), and (c): The Au leads are in the (001) 
direction, and there are 0, 2, and 1 additional Au atoms at the contacts, respectively, (d), (e), and (f): The Au leads are in 
(111) direction, and there are 0, 2, and 1 additional Au atoms at the contacts, respectively. The dashed line indicates the 
interface between the device region (C) and the left or right lead (L or R). 



Based on this method, here we report an investiga- 
tion of the molecular conductance of benzene connected 
to two Au leads with finite cross-section. We note that 
the effetcs of contact atomic structure in an unrelaxed S- 
anchored Au(lll) system have been investigated by Xue 
and Ratner using a cluster approach™^ We here focus 
on the effects of changing the atomic structure around 
the contacts, including the presence of an additional Au 
atom and changes in the molecule-lead separation, both 
of which are simple but important situations in break- 
junction experiments. In our calculation the effects of 
contact atomic relaxation and different lead orientations 
are fully considered. Also considered are different group- 
VI anchoring atoms. Our calculations show that the con- 
tact atomic structure is critically important to electron 
transport through the molecule: an additional Au atom 
at the contacts can increase the conductance by two or- 
ders of magnitude due to a LUMO-like resonance peak 
around the Fermi energy. We find that the dependence 
of the equilibrium conductance on the molecule-lead sep- 
aration can be complicated: either very weak or very 
strong depending on the separation regime. The reason- 
ability of the present systematic results are discussed in 
comparison with the results from Xue and RatnerSSl 



II. SYSTEMS INVESTIGATED AND 
COMPUTATIONAL DETAILS 

The systems we have studied consist of a benzene 
molecule connected to two Au leads of finite cross-section 
through a S atom located at the hollow site of the 



Au(001) or Au(lll) surface. The size of the Au leads in 
directions perpendicular to the transport direction is set 
to be 4^x4^ for the (001) lead and 3x3 for the (111) 
lead. To investigate the role of contact atomic structure, 
here we consider a very simple but possible situation: the 
presence of an additional Au atom at either one or both 
contacts, denoted by lAu and 2Au, respectively. We use 
the structural label (001)_lAu, for instance, to denote 
the system with the Au(001) leads and an additional Au 
atom at one of its contacts. To show the effect of change 
in the molecule- lead separation, we change rigidly (i.e., 
without any further structure relaxation) the contact Au- 
S distance (c?Au-s) m the vertical direction. The purpose 
of all these considerations is to simulate possible situ- 
ations in break-junction experiments in which different 
contact atomic structures may occur because of atomic 
fluctuations on the break surfaces and the molecule-lead 
separation can be adjusted by the MCB technique. To 
show the effect of contact atomic relaxation, we calculate 
the electron transmission for two independent cases: (1) 
Atoms in the leads are fixed at their bulk positions, and 
the molecule is fixed at the optimized structure of the 
isolated molecule with the dangling bond on the S atom 
saturated by an Au atom. The distance between the S 
atom and the Au surface, however, is optimized. This 
structure is called "unrelaxed" . (2) The structure of the 
molecule, the fist two atomic layers of the lead surfaces, 
as well as the molecule-lead separation are fully relaxed. 
The in-plane position of the S atom is, however, fixed at 
the hollow site. This structure is called "relaxed"; the 
structure of the relaxed systems are shown in Fig. 1. 

We use the efficient full DFT package SIESTA 41 to do 
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FIG. 2: Comparison between the transmission functions of the unrelaxed (dashed line) and the relaxed (solid line) systems. 
Different systems are indicated by their structure labels as defined in the text, (a), (b), (c), and (d) correspond to the atomic 
structures of (a), (b), (d), (e) in Fig. 1, respectively. 



the electronic structure calculation. It adopts a LCAO- 
like and finite-range numerical basis set and makes use of 
pseudopotentials for the atomic cores. We adopt a high 
level double zeta plus polarization (DZP) basis set for 
all atomic species. The PBE version of the generalized 
gradient approximation (GGA) 4 ^ is adopted for the elec- 
tron exchange and correlation, and optimized Troullier- 
Martins pseudopotentials 43 are used for the atomic cores. 
The atomic structure of the relaxed systems are opti- 
mized until the maximum residual force on all atoms is 
less than 0.02 eV/A. 

For the transport calculation, we divide an infinite 
LML system into three parts: left lead (L), right lead 
(i?), and device region (C) which contains the molecule 
and large parts of the left and right leads, as shown in 
Fig. 1, so that the molecule-lead interactions can be fully 
accommodated. Under a bias Vb the region C will be 
driven out of equilibrium. We have developed a simple 
while strict full self-consistent approach 27 to handle a 
steady state bias: The bias is included through the den- 
sity matrix of the region C (Dc) in the Green function 
calculation instead of the potential (He) in the DFT 
part. Specifically, we calculate Dc under the boundary 
condition that there is a potential difference V& between 
the left side of region C (together with the left lead) and 



the right side of C (together with the right lead), 



D C = -/ dE 



G C (E)T L (E + ^)Gl(E)f(E - Mi ) 



G C (E)T R (E - e -^)Gl(E)f(E - » R ) 



(1) 



where Gc(E) is the retarded Green function of region 
C (in which all the potential shifts are included) j2£ / 
is the Fermi function, and jiL and jiR the chemical po- 
tentials of the leads. Tl(E) and Tr(E) reflect the cou- 
pling at energy E between the C region and the leads L 
and R. The self-consistent loop is -> H C (DFT) -> G c 
Dc(NEGF) H c (DFT) ... until H c and D c 
converge™ The electron transmission through C is then 
related to Green functions by 



T(E,V h ) = Tr 



T L (E + e -^)G c (E)T R (E - ^)GUE) 

(2)" 

Note how Vb again appears in T here. Finally, the steady- 
state current is obtained by simply integrating the trans- 
mission over the energy window. 



eV b , 
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TABLE I: Calculated equilibrium conductance (G, in units 
of 2e 2 /h) and molecule- lead electron transfer (AQ, in units 
of electron, a positive value means that electrons are trans- 
ferred from lead to molecule) . Note the large effect of adding 
additional Au atoms. 
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III. RESULTS AND DISCUSSION 
A. Equilibrium Transmission and Electron Transfer 

In Fig. 2 we show the transmission functions for both 
the unrelaxed and the relaxed systems. The calculated 
values of equilibrium conductance are given in Table 1, 
together with the molecule-lead electron transfer (a pos- 
itive value means electrons are transfered from the Au 
leads to the molecule, including the two S atoms). 

Because we use the bulk Au structure for the leads in 
the unrelaxed cases, the contact atomic relaxation con- 
sists of two parts: (1) the relaxation of the bare Au 
lead with respect to its bulk structure, and (2) the relax- 
ation of both the leads and the molecule induced by the 
molecule-lead interaction. Our calculations of optimized 
atomic structure show that both parts are very small, as 
can be seen in Fig. 1. The small relaxation of the bare 
Au leads is consistent with the very small surface relax- 
ation of unreconstructed infinite Au surfaces. The small 
molecule-lead relaxation is understandable because the 
hollow site corresponds to a bulk atomic position, and 
so the absorption of a S atom will not markedly change 
the directional binding of the surface. Because the con- 
tact atomic relaxation is very small, its effect on electron 
transmission is only minor: as seen in Fig. 2 and Ta- 
ble 1, the induced change in equilibrium conductance is 
less than 100%. The very small molecule- lead relaxation 
justifies the reasonability of changing rigidly g?au-s for 
simulating the change in the molecule-lead separation. 

The results in Table 1 show that if there is no ad- 
ditional Au atom at the contacts (OAu systems), the 
Au(001) lead always yields a larger conductance (by 
about 6 times) than Au(lll). This can be understood 
by considering the contact atomic configuration in the 
two cases: the S atom touches four Au atoms on the 
Au(001) surface but only three on Au(lll). From Table 
1, the two OAu systems have very similar electron trans- 
fers (very small). In contrast, from Fig. 2 (a) and (c), 
their transmission functions around the Fermi energy are 
quite different: there is a small peak for Au(001) which is 
absent for Au(lll). Taken together, this indicates that 




FIG. 3: Local density of states within the energy window 
[-0.05, +0.05] eV around the Fermi energy for the systems 
(001)_0Au (upper) and (001)_2Au (lower). Note the channel 
with LUMO-like character in the latter which is absent in the 
former. 



the large difference in conductance between the two lead 
orientations is mainly due to the different molecule-lead 
couplings. Note that this large difference in conductance 
is significantly reduced by adding additional Au atoms 
to the contacts. This is obvious according to the above 
analysis: the additional Au atom reduces the structural 
difference for electron transport between the Au(001) and 
Au(lll) contacts. Another difference between the two 
lead orientations is that the overall structure of T(E) is 
totally different, as shown in Fig. 2: The T(E) functions 
of the Au(001) systems are quite smooth while those of 
the Au(lll) systems have many very sharp structures. 
We cannot understand well at this moment this differ- 
ence. Possible reasons include the too thin Au(lll) lead 
and its lower symmetry. This behavior of the Au(lll) 
lead has also been found in other calculations™ 

The introduction of Au atoms to the contacts causes a 
dramatic change in the conductance of the system. The 
calculated values of equilibrium conductance in Table 1 
show that an additional Au atom at both contacts in- 
creases the conductance by a factor of 14 for the Au(001) 
lead and 61 for Au(lll). From the calculated transmis- 
sion functions in Fig. 2, we see clearly that the two ad- 
ditional Au atoms produce a large resonance around the 
Fermi energy. It will be shown later that this originates 
from the LUMO level of the isolated molecule. The driv- 
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TABLE II: Calculated equilibrium conductance (G, in units 
of 2e 2 /h) and molecule- lead electron transfer ((AQ, in units 
of electron) for the Se- and Te- anchored systems. (The no- 
tations are the same as those in Table 1.) All the trends seen 
in the S anchored systems are also evident here. 
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ing forces causing the resonance peak are the molecule- 
lead electron transfer and coupling. As seen in Table 1, 
the introduction of the two additional Au atoms changes 
the sign of the electron transfer so that more electrons are 
transfered from the leads to the molecule, finally causing 
the (broadened) molecular LUMO level to line up with 
the chemical potential of the leads. This point will be 
confirmed directly later. To show more clearly the dif- 
ference between the systems with and without the addi- 
tional Au atoms we show in Fig. 3 the local density of 
states within the energy window [-0.05, +0.05] around the 
Fermi energy for the systems (001)_0Au and (001)_2Au. 
The two additional Au atoms clearly lead to a conduc- 
tance channel in the (001)_2Au system which is absent 
in the (001)_0Au case. Note that the spatial shape of the 
density states on the molecule indicates that the channel 
is formed from the LUMO level. When we add the Au to 
only one of the two contacts [see the atomic structures 
in Fig. 1 (c) and (f)], the increase in conductance is only 
minor (see the rows with label "lAu" in Table 1), espe- 
cially for the Au(001)-lead, indicating that the electron 
transmission through the molecule is then choked off by 
the less transparent contact. 



B. Other Anchoring Atoms 

In order to examine whether this behavior is common 
for other group- VI anchoring atoms, we carry out the 
same calculation for systems having similar geometry 
but anchored by Se and Te atoms. The results of equi- 
librium conductance and molecule-lead electron transfer 
are listed in Table 2. Clearly, the conclusions reached 
in the S-anchored systems also holds for the Se- and Te- 
anchored systems. 
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FIG. 4: Equilibrium conductance and lead-to-molecule elec- 
tron transfer as functions of the change in Au-S distance 
(AdAu-s), for (a) unrelaxed (001)_2Au system, (b) unrelaxed 
(001)_0Au system, and (c) relaxed (001)_lAu system. In (a) 
and (b) c?Au-s of the two contacts is rigidly and symmetri- 
cally changed, while in (c) only cZau-s of the contact without 
the additional Au is rigidly changed. The transmission func- 
tion for point C is shown in (d), while those for points A 
and B are already shown in Fig. 2(b) and (a), respectively. 
Note the large resonance conductance peak in panels (a)-(c) 
accompanied by significant electron transfer to the molecule. 
The similarity of (d) to Fig. 2(b) shows that increased Au-S 
separation has an effect comparable to that of the extra Au 
atoms. 



C. Dependence on Molecule-Lead Separation 

In break-junction experiments the molecule-lead sepa- 
ration may not be at its equilibrium value but rather is 
probably lengthened or compressed because of the mis- 
match between the molecular length and the junction 
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break. To simulate this situation, here we calculate the 
equilibrium conductance as a function of the change in 
<^Au-s (A^Au-s) with regard to its equilibrium value for 
three systems: unrelaxed (001)_2Au, (001)_0Au, and re- 
laxed (001)_lAu. For the first two cases, c?Au-s of both 
contacts will be changed rigidly while maintaining the 
symmetry of the system, while for the third case only 
<^Au-s of the contact without the additional Au atom 
will be changed rigidly. As mentioned previously, the 
very small molecule- lead relaxation justifies this treat- 
ment. The results are shown in Fig.0] 

There is a large resonance peak in the conductance 
curve for all three systems. For the (001)_2Au case, the 
equilibrium Au-S distance is very close to the position of 
the resonance peak, while for the other two systems the 
equilibrium Au-S distance is about 1.4A away from the 
position of the resonance peak. Along with the increase 
of Ag^au-s the amount of electron transferred from the 
leads to the molecule increases and reach its maximum 
around the resonance peak. If we assume the mechanism 
producing the resonance is the same for the three sys- 
tems, then T(E) at point C in Fig.^b) should be similar 
to that at point A in Fig. Ufa), even though transmission 
at points A and B are very different [see Fig. 2(b) and 
(a)]. To check our point we plot the transmission func- 
tion for point C in Fig.0Jd). It is clear that, as expected, 
this T(E) is quite similar to that of Fig. 2(b). 

Before continuing to investigate the mechanism of the 
resonance, we briefly pause to compare to a previous cal- 
culation: the S-anchored Au(lll)-lead system was inves- 
tigated using a cluster methociS2i2& which did not, how- 
ever, include the atomic relaxation of the contact. In the 
area where they overlap, there is qualitative agreement 
between the two calculations; that is, they both find that 
the presence of the additional Au atoms at the contacts 
significantly increases the equilibrium conductance and 
that this effect is quite similar to that from increasing 
the contact Au-S distance. However, there are significant 
quantitative differences between the two sets of results. 
For instance, the gap in T(E) for the (lll)_0Au case in 
Ref. 29 is much larger than our result [compare Fig. 9 
in Ref. |2£j to Fig. 2(c) here]; the latter is comparable to 
the GGA HOMO-LUMO gap of the molecule (~ 0.3eV) 
and thus appear more reasonable physically. These sig- 
nificant differences are probably caused by the very dif- 
ferent techniques used in the two calculations: (1) In our 
method the electronic structure of the molecule and the 
leads are treated exactly on the same footingdA while in 
Ref.HiU the molecule is treated by DFT but the leads 
are treated by a tight-binding approach. (2) We use pe- 
riodic boundary conditions with large parts of the leads 
included in the device region; as a result, the semi-infinite 
leads and the infinite LML system are nearly perfect in 
geometry^ In contrast, in Ref. 2903(1 a cluster geometry 
is used for the device region, introducing some surface ef- 
fects. Our consistent treatment of the semi-infinite leads 
also allows us to analyze further the mechanism of the 
resonance peaks, as shown in the next section. 
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FIG. 5: Contour plots of transmission as a function of both 
energy and Ao?Au-s for (a) unrelaxed (001)_2Au system, (b) 
unrelaxed (001)_0Au system, and (c) relaxed (001)_lAu sys- 
tem. The meaning of Ag?au-s is the same as in Fig. 0] The 
contributions from the HOMO and LUMO orbitals of the iso- 
lated molecule are indicated in (a) and (b). 



D. Resonance Mechanism 

In order to show more clearly the mechanism of the 
resonance peak around the Fermi energy and the mod- 
ification in transport properties caused by changing the 
contact Au-S distance, we show in Fig. 5 contour plots 
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of transmission coefficient as functions of energy and 
A^Au-s f° r the three S-anchored systems, (001)_2Au, 
(001)_0Au, and (001)_lAu. For the (001)_2Au system, 
when the Au-S distance is large (A^Au-s > 1-OA) there 
are three energies in the energy window which contribute 
to the electron transmission. Comparison with the level 
structure of the isolated molecule shows that these three 
energies correspond to the LUMO, HOMO, and HOMO- 
1 of the isolated molecule, respectively. As AoIau-s de- 
creases, the increasing molecule-lead coupling broadens 
these three levels. The two HOMO states shift gradually 
to lower energies as well. It is clear that the resonance 
peak around A^Au-s = 0.2A in Fig.^a) originates from 
the (broadened) LUMO contribution. 

For the (001)_0Au case and large Au-S distance 
(A^Au-s > 1-4A), the result is completely similar to 
that for the (001)_2Au system. This is a further indi- 
cation that the two resonance peaks in Fig. 0] (a) and 
(b) have the same character, as we already argued from 
the similarity similarity T(E) in Fig. 2(b) and Fig.^d). 
The role of the two additional Au atoms at the contacts 
is equivalent to that from increasing the surface-S dis- 
tance. As A<iAu-s is decreased to smaller than 1.0 A the 
strong molecule-lead coupling changes significantly the 
local electronic structure, and we cannot distinguish the 
individual contribution from the molecular orbitals any- 
more. Finally, at A^Au-s = O.OA the transmission func- 
tion becomes totally different from that for large At/Au-s, 
as shown in Fig. 2(a). 

An interesting thing we should notice in Fig. 5(b) is 
that within the large range of A(iAu— s ^ 

0-lA the equi- 
librium conductance is actually very insensitive to the 
Au-S distance. This indicates that the molecule-lead sep- 
aration dependence of the equilibrium conductance can 
be quite complicated: it can be either very strong (for 
large A^Au-s) or very weak (for small A^Au-s)- 

For the (001)_lAu system, because of the strong cou- 
pling on the left side [see Fig. 1 (c)] we cannot recognize 
the individual contributions from the molecular orbitals 
in Fig. 5(b) even for large A^Au-s- However, there is a 
similar LUMO-like contribution to the resonance peak in 
Fig.^c), and the equilibrium conductance is also insen- 
sitive to the Au-S distance for Ag£au-s ~ 0-lA. This in- 
dicates that the total electron transmission is dominated 
by the weakly coupled contact. 



E. I-V Curve 

The large resonance peak in the equilibrium transmis- 
sion function around the Fermi energy implies a large 
negative differential conductance when a bias voltage is 
applied. We would like to show this explicitly. In order 
to avoid the large computational effort for the I—V curve 
of this system, we with a similar but slightly smaller sys- 
tem: two smaller 2^/2 x 2\/2 Au(001) leads instead. The 
structure and equilibrium transmission function of the 
small system are shown in Fig. 6. It can been seen that 
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Bias (V) 

FIG. 6: I-V curve of the smaller (001)_2Au system shown in 
the upper inset, whose transmission function under zero bias 
is shown in the lower inset. Note that there is a large reso- 
nance peak in T(E) around the Fermi energy which causes 
a large negative differential conductance in the I-V curve 
around V b = 0.2-0.5 V. 



its transmission function is somewhat different from that 
of the large system, but the feature of the large resonance 
peak around the Fermi energy is the same. The calcu- 
lated I — V curve given in Fig. 6 shows clearly a large 
negative differential conductance around Vb ~ 0.2 - 0.5 
V. 



IV. SUMMARY 

By using a density functional theory calculation for 
molecular electronic structure and a Green function 
method for electron transport, we have calculated from 
first principles the molecular conductance of benzene 
sandwiched between two Au leads in different ways. In 
our calculation, the effects of contact atomic relaxation, 
two different lead orientations, and different anchoring 
atoms are fully considered. We focused on the effects 
of the change in atomic structure around the contacts, 
including the presence of an additional Au atom and 
changes in the molecule- lead separation, as an effort 
to simulate possible situations in break-junction exper- 
iments. Our calculations clearly show the critical role 
of the contact atomic structure in electronic transport 
through molecules: 

(1) The presence of an additional Au atom at each 
of the two contacts can increase the equilibrium conduc- 
tance by one to two orders of magnitude regardless of the 
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contact atomic structure or group- VI anchoring atom. 
The mechanism is the creation of a LUMO-like resonance 
peak around the Fermi energy, which also leads to nega- 
tive differential conductance under applied bias. 

(2) The presence of the additional Au atom at only one 
contact will give only a minor increase in conductance 
because the electron transmission is then choked off by 
the other non-transparent contact. 

(3) Because of the different molecule-lead coupling, the 
Au(OOl) lead will lead to an equilibrium conductance up 
to 6 times larger than that of the Au(lll) lead. This dif- 
ference is significantly reduced by adding the additional 



Au atom at the contacts. 

(4) The dependence of the equilibrium conductance on 
the molecule- lead separation is subtle: it can range from 
either very weak to very strong depending on the sepa- 
ration regimes. 
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